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Abstract 

The measurement of the average depth of the shower maximum is the 

most commonly used observable for the possible inference of the primary 

i-C \ cosmic-ray mass composition. Currently, different experimental Collabora- 

I I tions process and present their data not in the same way, leading to problems 

2 ■ in the comparability and interpretation of the results. Whereas (Xmax) is ex- 

'^ • pected to be proportional to (In A) in ideal conditions, we demonstrate that 

•^ the finite field-of-view of fluorescence telescopes plus the attenuation in the 

atmosphere introduce a non-linearity into this relation, which is specific for 

^ \ each particular detector setup. 

Keywords: cosmic-rays, extensive air-showers, shower maximum, telescope. 



OO . 

C"^ '. field-of-view 



m 

O ' 1. Introduction 

(N ; 

Ultra-high energy cosmic-ray (UHECR) particles are the highest energy 

. !^ ■ particles observed by humankind, well beyond what is accessible at the Large 

/\ . Hadron Collider. Current experiments [l|, Ij, |3| are recording large amounts 

o3 ! of high quality data. Recent observations do not draw a simple nor consistent 

picture of the nature of UHECR particles. There are hints of anisotropy in the 

arrival directions [^ , favouring the presence of light primary particles. Using 

the current high energy hadronic interaction models, typical air-shower ob- 

servables do not generally support the light component hypothesis: the Pierre 
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Auger Collaboration observes the depth of the electromagnetic shower maxi- 
mum, muon production maximum, rise-time of the shower front, etc. that are 
all pointing in the direction of a primary composition dominated by heavy 
particles at ultra-high energies J5|, |6|. On the other hand, the HiRes Collab- 
oration claims that their measurement of the depth of the shower maximum 
is compatible with protons up to the highest energies |7|, which is currently 
supported by the first TA data [sl, |9| . It is important to notice that whereas 
Auger publishes (Xmax) in the atmosphere by using a fiducial volume selec- 
tion to minimize acceptance biases, HiRes and TA presently do not apply 
such corrections and present (Xmax) at the detector level. Also the possible 



difference of the energy scale between Auger and HiRes/TA (^25 %) [10 



contribute to the difficulties to explain the data [ll|, Il2[ . 

The current situation could point to exciting physics aspects that are 
about to emerge with higher accumulated event statistics, as for instance, 
systematic differences of the cosmic-ray flux on the northern and southern 
hemisphere. The Pierre Auger Observatory observes the southern, while 
HiRes/TA the norther sky, with very different astrophysical objects in their 
direct field-of-view. It is clear, that the understanding of the primary mass 
composition of UHECR is one of the major steps towards the final solution 
of the UHECR puzzle. Depending on a reliable measurement of the mass 
composition very different scenarios of the nature of UHECR will finally 
emerge. 

In this article we discuss the concept of observed versus true Xmax- 
distributions to finally demonstrate that the observed average of depth of 
shower maxima, (Xmax), can result in a non-linear relation with the average 
logarithm of the mass number, (In A) , of the UHECR. These effects are very 
specific for a given experimental setup and have to be accounted for in order 
to allow a comparison or interpretation of these data. 

2. Extensive Air-Showers and Fluorescence Telescope 

Fluorescence telescopes measure the ultraviolet light, which is propor- 
tional to the energy deposited by the passage of the charged particles of 
the air-shower cascade through the atmosphere. This makes possible to re- 



construct the longitudinal profile of the electromagnetic shower [13[. The 
atmospheric depth at which the energy deposit is maximal is the shower 
maximum, X^ax, and is related to the nature of the primary particle. At the 
same primary energy, primary nuclei with mass A produce air-showers with 



shower maxima at different atmospheric depths. This can be approximated 



as 
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771 
-^max ~ Aint + In 2 Xq In — , 

where Ajnt is the cross-section of cosmic-ray primaries with air and Xq ~ 
37g/cm^ is the electromagnetic radiation length in air. The structure of this 



equation, X^ax = a+6 In E/A, holds also for more general considerations [15 
It has been typically used to compute the average logarithmic mass from the 
data of the average depth of the shower maximum 

{\nA) = -[a + blnE- (X^^x)] = ^ [(X^ax)p - (X^ax)] • 

Thus, there is a linear dependence of the measured average shower maximum 
from the average mass of the cosmic-ray primary particles. 

Instead of calculating (In A) one can also compute the component frac- 
tions in a model with two primary cosmic-ray species. This is typically done 
for the proton fraction / under the assumption of a simple proton/iron mix- 
ture. For this case the relation of (Xmax) is 

(-^max) = /(-^max/p + (1 ~ /) (-^max/Fe, (1) 

which is linear in /. 

Thus, the distinction between an admixture of different primaries or a 
pure primary of an intermediate mass, are indistinguishable on the base of a 
(-^max) value alone. The resulting equivalence is 

(In A) = ^ [(X^ax)p - (X^ax)Fc] (1 - /)• 

3. The effect of the telescope acceptance 

The objective of this paper is not to reproduce any particular detec- 
tor configuration, but to demonstrate that the relation between (X^ax) and 
(In A) depends critically on the detector acceptance, which varies with the 
primary energy. Thus, the telescope acceptance has an impact on the in- 
terpretation of (Xmax) data, and in general any other momentum of the 
Xmax-distribution (see also |l6|). 
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Figure 1: Left panel: Geometry of our simulation study. Right panel: Resulting efficiency, 
^obs/-/Vgon, depending on the distance of air-shower cores from the telescope. 



Fluorescence telescopes have a limited viewing angle defined by the optics. 
Showers falling outside the field-of-view cannot be detected. But even for 
observed showers the shower maximum can be outside the field-of-view of 
the telescope and can then not be reconstructed reliably. The angular field- 
of-view defines a volume in the atmosphere where showers can effectively be 
fully reconstructed. One of the most important criteria of this is that the 
shower maximum is located within the field-of-view. This volume is depicted 
in Fig. [1] (left). When showers fall very close to the telescope, the transverse 
area of the field-of-view is small, meaning that many shallow and/or deep 
showers are not fully reconstructed. This effect is naturally most severe 
in the energy range close to the lower detection threshold for air-showers 
where the distance to the events is limited by the small amount of generated 
fluorescence light. For air-showers detected at larger distance to the telescope 
the geometrical field-of-view cone is much larger. However, fluorescence light 
can travel only a limited distance in the atmosphere before being absorbed. 
At some point, light emitted at the shower axis is attenuated too much 
and cannot be observed any more. This is responsible for the characteristic 
rapid drop of the efficiency as shown in Fig. [1] (right) beyond 20 km. The 
Rayleigh absorption rate is proportional to the density of the atmosphere, 
and thus depends on changing air pressure as well as on the height in the 
atmosphere at which a particular shower is developing. In this work we do 
not discuss aerosol related absorption. It is much smaller than Rayleigh 
absorption and to include it here does not help the argument, but just adds 
additional complexity. 



By means of a toy Monte Carlo we reproduce all parameters relevant 
for the modelling of air-shower observation with fluorescence telescopes at 
an arbitrary primary energy. Our simulations focus on the geometric and 
atmospheric effects and simplifies the situation as much as possible. We do 
not attempt to reproduce a specific detector setup, but only want to show 
the inevitable impact on the interpretation of the observed data. 

The geometry of our simulation is shown in Fig. [T] The vertical depth of 
the telescope is at ~ 880g/cm^, which corresponds to h=1.5km, as explained 
below, and the opening angle of the field-of-view is 20°. We consider air- 
showers of fixed primary energy, with random arrival directions, sampled 
from dA^/dcos6' =const up to 60° zenith angle, as well as distances / from 
the telescope, sampled from dN/dl oc /. The requirement of a particular 
shower to be observed is that at least A^det = 100 photons reach the telescope 
from the location of the shower maximum on the shower axis. The number 
of photons is computed as 

iVdet = iVph^dia/r' exp ( - -^ ) , (2) 

where A^ph is the number of fluorescence photons emitted at the shower max- 
imum, Adia = 10 m^ is the aperture of the telescope, r is the geometric dis- 
tance, t = j p{z)dr the integrated depth distance from the telescope to the 
location of ^Y'max, and Aabs = 1000 g/cm^ is the photon absorption length. For 
air showers of primary energy Eq and a fluorescence yield of Ifluo = 5/MeV 
the number of photons emitted at X^ax is 

iVph=(^^j Ifluo pc At, (3) 

where At = 100 ns is the telescope sampling time and (di?/dX)max is the 
energy deposit at the shower maximum, which is lO~^^-'^'^+^s(Eo/eV) GeVcm^/g 
in very good approximatioiij. For air showers with Eq = lO^*^'^ eV this yields 
A^ph ~ 2 ■ 10^*^, which is the default value throughout this paper if not stated 
otherwise. The atmospheric density profile p{z) is exponential with a scale 
height of 8 km and a pressure at the height above sea level, h = 1.5 km, of 
the telescope of ~ 880g/cm^. 



This is obtained with the SIBYLL [17| interaction model. 
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Figure 2: Efficiency e(Xniax) of tlie X,„ax reconstruction as a function of Xmax for a 
particular detector, as calculated by our toy Monte Carlo. 



When we generate events with a flat distribution of Xmax, we find that air- 
showers with different Xmax are observed with different efficiency e(Xmax) = 
Nohs/Ngen- In Fig. |2]we show this efficiency as a function of X^ax- The peak 
of this distribution is related to the average atmospheric slant depth in the 
volume observed by the telescope. Every telescope setup has limits both for 
very small as well as very large values of X^^x (see e.g. [l8|). While large 
integrated atmospheric depths can in principle always be achieved by very 
inclined geometries, there is a strict bound on the minimal observed depth 
even for vertical events, which is related to the maximum possible observation 
distance and the upper elevation boundary of the field-of-view. This effect 
becomes more relevant for lower energy air-showers, since here the maximum 
observation distance is smaller. 

Given the average efficiency e(Xmax) of a telescope setup, the measured 
distribution of Xmax is related to the parent distribution via 

dN \ ,^ . f AN 



AY I ^ ' \ r\ Y 



Small changes in the telescope setup can yield very different e(Xinax)- If also 
energy, Xmax resolution, and possible Xmax reconstruction biases effects were 
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Figure 3: The two input distributions corresponding to light and heavy primary cosmic- 
ray particles in our study. The distributions are Exponential convoluted with a Gaussian. 
The used parameters of Eq. (|4]) are /i = 815g/cm^, a = 50g/cm^ and r = 45g/cm^ for 
the light component as well b& p. — 1005 g/cm^, ct = 30g/cm^ and r = 5g/cm^ for the 
heavy component. 



included the relation between the measured and parent distribution would 
become more complex 19|, l20|, but this does not help the clarity of our 
argument. 

In the following we demonstrate how the telescope acceptance, e(Xmax), 
affects in particular also simple analyses as for example related to (Xyaax)- 
For this purpose we generate input X^ax distributions derived from an Ex- 
ponential convoluted with a Gaussian function 



dN 

dXma, 



A^, 



cvt 



2r 



ga2/(2r2) g(/.-X^ax 



)A 



erfc 



Ai - -^max - cr'^/r 



V^i 



a 



(4) 



In total this function has three shape parameters: the mean and width of 
the Gaussian, fi and a, and the exponential slope, r. We generate two 
distributions (c.f. Fig. [3]), the first one corresponding to light cosmic-ray 
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Figure 4: Dependence of the bias A(X,„ax) from the energy and thus the distance of air- 
showers from the telescope. The upper axis indicates the resulting average distance of 
showers from the telescope depending on the primary energy. 



primaries has an average of 750 g/cm^ and an RMS of 66.5 g/cm^, the second 
one corresponds to heavy cosmic-ray primaries has an average of 650 g/cm^ 
and an RMS of 30g/cm^. These values are chosen since they correspond 
roughly to typical values for proton or iron induced air-showers respectively. 
The average efficiencies 
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diV 
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e(Xjnax)dXn 



for the observation of X^ax with our telescope setup for these two distri- 
butions are eught = 33.1% and eheavy = 16.7%. The average X^ax ob- 
served by the telescope setup are (Xinax)iight — 770g/cm^ for the light and 
(-^max)heavy = 660 g/cm^ for the heavy component, which is 10 — 20g/cm^ 
biased with respect to the input distributions. The bias A(Xmax), de- 
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Figure 5: The average X^ax as a function of the hght cosmic-ray fraction, /. The observed 
(^max) values for pure hght and pure heavy compositions are shifted with respect to the 
true (Xinax)- The relation between (Xmax) and /, and therefore with (In A) by Eq. [l] is 
linear for the true v\",nax distributions, whereas it is non-linear for the case of the observed 
(-''^max)- The curve labeled "Observation A" indicates the effect on a pure composition of 
intermediate mass. 



fined as the difference between the observed and true (Xmax) value, changes 
with distance to the telescope and thus also with the primary energy of the 
cosmic-ray particles according to Eqs. ([2]) and (|3]). In Fig. Hlthis is shown for 
the cases of pure light, pure heavy and an equal mixture of light and heavy 
primaries. For this study we consider a change of (Xmax) with energy of 
d(Xmax)/dlg-E = 60g/cm^ with respect to the default values used otherwise 
in this paper. It is interesting to note that the resulting bias depends on the 
underlying mass composition. In general, for smaller energies the showers are 
located closer to the telescope and the effect becomes stronger. For exam- 
ple, the low energy fluorescence telescope extensions HEAT 21| at the Pierre 
Auger Observatory and TALE 221] at the Telescope Array are important to 
limit these biases by providing a much wider fleld-of-view for showers at close 
distances. Furthermore, the bias can be positive as well as negative, depend- 
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Figure 6: Impact of different values of r = eiight/chcavy on the non-linearity. 



ing on the overlap of the true X^ax distribution with the telescope efficiency 
function, c.f. Fig. |2l However, the impact will qualitatively be always as 
shown in Fig. HI at short distances the field-of-view cuts shallow showers and 
thus (Xmax) is Overestimated, at larger distances this effects becomes smaller 
and eventually might even reverse leading to an underestimation of the true 
(^max) value. It is important to realize that the primary effect is due to a 
different distance of showers, and that any parameter that affects the typical 
observation distance will have a similar impact. Such parameters are typi- 
cally related to the detector setup, for example the diaphragm opening Ajja, 
the sampling time At, but also the atmospheric density profile and optical 
absorption characteristics. 

If we consider a mixture of light and heavy particles at fixed primary 
energy with a given fraction / of light and 1 — / of heavy primary particles, 
the average observed shower maximum becomes 



X^rnax/ 



/'"(-'^ max /light + (1 /)\-^ max) heavy 



(5) 



where r 



fr + l-f 
ehght/^heavy This relation is non-linear in / and thus not in In A 
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In Figure we show the results of this study. The true as well as the observed 
(-^max) are shown, together with the linear interpolation from Eq. ([T]) and 
the non-linear version Eq. (|5]). The latter describes perfectly well the results 
of our Monte Carlo study. 

We also study how a different detector setup affects our result. In general, 
by changing parameters of the telescope detector setup, very different values 
of r can be obtained. In Fig. Elwe show how the predictions of Eq. i^ change 
for a wide range of different values of r. It summarizes the effect described 
in this paper, and the miss-interpretation it can induce on the inference of 
the primary mass composition from (Xmax)- For instance, r = 3 would imply 
that a real fraction 50 % proton and 50 % iron, that is f=0.5 and (In A) = 2, 
would be interpreted as 75% proton and 25% iron, {In A) = 1, if just the 
linear relation Eq. ([1]) is assumed. 

4. Summary 

The data collected by a fluorescence telescope are affected by acceptance 
effects, which can have an important impact on the detailed interpretation 
of the data. This is problematic for researchers not from within a particular 
experimental collaboration who, for example, have no access to a detailed 
Monte Carlo simulation of the detector setup. 

The average efficiency e(Xmax), which describes the response to a fiat dis- 
tribution in Xjnax, is one of the crucial properties of a fluorescence telescope. 
However, in addition to the acceptance effects also the detector resolution 
and Xmax reconstruction biases have to be known for a full data analysis. 

We have demonstrated that the bias in the average shower maximum in- 
troduced by the acceptance of fluorescence telescopes induces a non linearity 
between (X^ax) and {In A). 
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